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Abstract 

Isotopic fluctuations in fragment formation are investigated in a quasi-analytical description 
of the spinodal decomposition scenario. By exploiting the fluctuation-dissipation relations the 
covariance matrix of density fluctuations is derived as a function of the wave vector k for nuclear 
matter at given values of density, charge asymmetry, temperature, and of the time that the system 
spends in the instability region. Then density fluctuations in ordinary space are implemented with 
a Fourier transform performed in a flnite cubic lattice. Inside this box, domains with different 
density coexist, from which clusters of nucleons eventually emerge. Within our approach, the 
isotopic distributions are determined by the N/Z ratio of the leading unstable isoscalarlike mode 
and by isovectorlike fluctuations present in the matter undergoing the spinodal decomposition. 
Hence the average value of the N/Z ratio of clusters and the width of the relative distribution reflect 
the properties of the symmetry energy. Generating a large number of events, these calculations 
allow a careful investigation of the cluster isotopic content as a function of the cluster density. A 
uniform decrease of the average charge asymmetry and of the width of the isotopic distributions 
with increasing density is observed. Finally we remark that the results essentially refer to the early 
break-up of the system. 
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I. INTRODUCTION 



Reactions with charge-asymmetric systems open the possibihty to learn about the proper- 
ties of the symmetry term of the nuclear interaction in conditions of density and temperature 
away from the ordinary values. In particular, the study of multi-fragmentation mechanisms 
in neutron-rich systems should allow to get information on the behaviour of the symmetry 
energy in density regions below saturation, where the nuclear system may undergo a liquid- 
gas phase transition. Constraints on the form of the density dependence of the symmetry 
energy are important not only for a better knowledge of the nucleon — nucleon interaction, 
and hence its extrapolation to the structure of exotic nuclei [l|, 0, 0] , but also for the study 
of the neutron star crust, and of supernova explosions, where a key issue is the clustering of 
low-density matter [1, H, [g], 0] . 

The dynamics of first-order phase transitions is often induced by instabilities against 
fluctuations of the order parameter. In dissipative heavy ion collisions, nuclear matter may 
be pushed inside the coexistence region of the nuclear liquid-gas phase diagram. Then, 
the observed abundant fragment formation may take place through a rapid amplification of 
spinodal instabilities. Experimental results pleading in favor of such a spinodal decomposi- 
tion have recently been reported jl, 0] . Spinodal instabilities in charge-asymmetric systems 
have been widely investigated from the theoretical point of view. The most important effect 
induced by the charge asymmetry is the so-called isospin distillation: fragments (liquid) ap- 
pear more symmetric with respect to the initial matter, while light particles (gas) are more 
neutron-rich [l^, 11, 12, 13, 14, 15|. The amplitude of this effect depends on specific prop- 
erties of the isovector part of the nuclear interaction, namely on the value and the derivative 
of the symmetry energy at low density [l6|. Moreover, apart from the distillation effect, that 
determines the average fragment isotopic composition, the symmetry energy value influences 
also the width of the isotopic distributions. This feature has been recently exploited in the 
so-called isoscaling analysis [13], where information on the symmetry energy behaviour is 
extracted from the study of ratio of the isotope yields obtained from two reactions with 



different charge asymmetry [18|, Il9|, |20|, [21 



In this paper we focus on a detailed study of isotopic properties of nucleon clusters, as 
obtained within the spinodal decomposition scenario. In order to select this mechanism, we 
consider nuclear matter initialized at a given temperature and at low density, inside a box 
with periodic boundary conditions, under the action of a stochastic field self-consistently 
determined [12, 2^. In this case, the normal modes of the density fluctuations are plane 



waves. The instabilities are treated in the linear approximation, i.e. retaining only first- 
order terms in the stochastic field or in the density fluctuations. This leads to a quasi- 
analytical description of the growth of instabilities and the consequent formation of clusters 
of nucleons, that can be considered as excited primary fragments. For each considered 
system, a large variety of possible outcomes may be obtained, according to the initial density 
fluctuation values. In this way it is possible to collect a large number of events with much 
reduced computational effort, allowing to perform a thorough analysis of the isotopic content 
of nucleon clusters, the isospin distillation and isotopic distributions, in connection with the 
ingredients of the effective nuclear interaction employed. 

We would like to stress that our analysis refers essentially to the early break-up of the 
system, where excited nucleon clusters, reflecting the nuclear matter properties, can be 
recognized. At this level, shell effects, that influence the N/Z ratio of real nuclei, are not 
considered. In particular, we will discuss the relation between isotopic properties and density 
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of the primary fragments (at the moment of fragment formation). In fact, in the spinodal 
decomposition scenario, domains with different density values, from which clusters of nu- 
cleons eventually emerge, coexist. Clusters of intermediate mass, that may originate from 
the vapour phase or from the liquid phase as well, may have different isotopic properties, 
depending on the density value of the domain where they are formed. For instance, the 
isospin distillation mechanism, that is a feature of isoscalarlike oscillations, is more and 
more effective as soon as the density gets larger in the domain considered. 

In the actual disassembly of a nuclear system (in central nuclear collisions) the vapour 
phase is generally formed at the surface of the system, where also a larger radial collective 
flow is observed. Hence one could expect a relation between fragment isotopic observables 
and kinematical properties [23]. In particular, as we will discuss in the following, clusters 
emerging from low density regions should be more neutron-rich and with wider isotopic 
distributions. These features are peculiar of the spinodal decomposition scenario and of the 
occurrence of first-order liquid-gas transitions. Hence, in addition to the information that 
can be gained about the low-density behaviour of the symmetry energy, this investigation 
would allow to shed some light on the fragmentation mechanism itself. 

The paper is organized as it follows. In Sec. II we outline the formalism developed in 



Refs. |22l . |23| for infinite nuclear matter and its implementation in a finite cubic lattice. 



In Sec. HI we discuss the results of our calculations. Finally, in Sec. IV a summary and 
conclusions are given. 



II. FORMALISM 

In this section we outline the main steps of the general formalism developed in Refs. 



[22l . |23j to evaluate fluctuations of the one-body density for asymmetric nuclear matter 
inside the unstable (spinodal) region of the nuclear matter phase diagram. In our approach 
fluctuations of the proton and neutron densities are induced by a stochastic field which 
couples with the costituents of the system. By means of the fluctuation-dissipation theorem 
the stochastic field is self-consistently determined. For nuclear matter in the spinodal region 
the time growth of fluctuations is essentially due to the unstable mean field. Therefore 
we neglect the effects of nucleon-nucleon collisions in the time evolution of fluctuations. 
Collisions would mainly add a damping to the growth rate of the fluctuations and should 
not change the main results of our calculations, at least at a qualitative level. 

The approach of Refs. 22, 2^ has reference to infinite nuclear matter, thus normal 
modes are plane waves associated with wave numbers k. The growth of fluctuations leads 
to the formation of large density domains, that can be associated with fragments. In Refs. 
[2^ . 23| the fragment recognition is based on a procedure which allows to determine the 



probability distribution of the domains containing correlated density fluctuations. We have 
then identified the pattern of correlated domains with the fragmentation pattern. In such 
a way, we could make predictions on the distributions of clusters. This procedure is based 
on an "ansatz" which relates the Gaussian distributions for the different modes of density 
fluctuations to the probability distribution of the correlated domains in the ordinary space. 

In the present paper we follow a different strategy. From the probability distributions of 
fluctuations in the k space, induced by the stochastic field, we directly generate a certain 
number of events characterized by stochastic distributions of density fluctuations in the 
coordinate space. This is accomplished by performing a Fourier expansion in a finite cubic 
lattice of volume V. 
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In this way a substantial improvement of the approach of Refs. 22|, |23|] is obtained. 
Event by event analyses can be performed. Moreover, one can investigate new features of 
the fragmentation process, such as, for instance, the isotopic content of nucleon clusters as 
a function of their density. 



A. Density fluctuations 

The key quantity for the evaluation of the density fluctuations is the density-density 
response function. In a linear approximation for the stochastic field, the Fourier transform 



of the response function is given by the equation [22, [23 



Di^,{k,u) = Df\k,uj)5i^j + T.iDf\k,uj)Au{k)Di^,{k,uj) (1) 

where Df'\k,uj) is the non-interacting particle-hole propagator and Ai^i{k) is the Fourier 
transform of the nucleon-nucleon effective interaction. The subscripts of the various quan- 
tities take the values 1 and 2 for protons and neutrons respectively. In this paper we use 
units such that h = c = ks = 1- 

In asymmetric nuclear matter isovector and isoscalar fluctuations are coupled. However 
one can still separate oscillations with neutrons and protons moving in phase (isoscalarlike) 
or out of phase (isovectorlike) and add the contributions of the corresponding variances. 
This can be done because the time scales of isoscalarlike oscillations and of isovectorlike 
oscillations are very different for nuclear matter with values of temperature and density 
sufficiently close to the borders of the spinodal region. Indeed, the growth rate of the 
unstable isoscalarlike modes, F^, turns out to be much smaller than the real frequencies, 
of isovectorlike modes 22, 2^. Thus the two kinds of fluctuations can be considered 



two independent stochastic processes. In addition, for the relevant values of the magnitude 
of the wave vector k, the following inequalities F^/T < 1 and /T > 1 hold. 

As a consequence of the linear approximation for the stochastic field, the probability 
distribution of density fluctuations, P[Sgi(k,t)], is given by a product of Gaussian distribu- 
tions: 

P[<5ft(k,t)] = iVexp ( - i-^^5f,^(k*,t)(a2(A;,t))r;5^^.(k,t)) , (2) 

k i,j 

where is a normalization constant. Each single factor corresponds to a stochastic process 



for a given wave vector k |25| . The covariance matrix for the isoscalarlike fluctuations 

aUk,t)=TC-{k)^^{e'^^'-l), (3) 

can be evaluated taking the classical (thermodynamical) limit uj/T <^ 1 of the fluctuation- 
dissipation relation. The coefficients Cl'^j{k) are the residues, times i, of the response function 
at the pole u = iT^ j22|. 

Since the time scale of the isovectorlike fluctuations is much shorter than that of the 
growth of the unstable modes, for the covariance matrix of the isovectorlike fluctuations we 
can take its asymptotic value for t cxd, with nuclear matter at given values of density, 
temperature and asymmetry, 

^U^ = Ct:^ik), (4) 
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where Clj{k) represent the residues of the response function at the real pole u^jy. This 
equation has been obtained by exploiting the fluctuation-dissipation theorem, now in the 



limit u/T >> 1 [22 



It should be noticed that, in the low temperature limit, the isovector fluctuation amplitude 
is essentially different from the value expected at the thermodynamical limit, u << T, where 
the fluctuation variance is proportional to the temperature and inversely related to the value 
of the symmetry energy [26| . 

With the aim to preserve a simple formalism, which, to some extent, allows us to per- 
form calculations analytically, we have evaluated the response function of Eq. ([1]) within 
a semiclassical approximation (Vlasov equation). However, this approximation cannot be 
extended to values oi k > kp [kp being the neutron or proton Fermi momentum). We 
cure this shortcoming with an ad hoc specific recipe. We calculate poles and residues of the 
response function within the Vlasov approximation for k ^ 1.1 /cjvf? where fcjv/ is defined by 
^k>kM < 0. With the physical parameters used in this work /cm is about the neutron Fermi 
momentum. For k > 1.1 /c^ we take for the total variance of density fluctuation distributions 



its value for /c — > oo, af j{k,t) = SijQt |27|, with Qi being the density of ^-species nucleons. 
This recipe is suggested by two circumstances. Firstly, for kM < k < 1.1 kM the variance, 
which relaxes towards its asymptotic value since Ffc < 0, remains appreciably larger than 
its limit Qi for k oo. Secondly, the latter is approached already for values of k slightly 
larger than I.I/cm, as it can be shown by explicit calculations including quantum effects. 
The non-accurate evaluation of the response function in the interval 1.1 kM < k ^ 1.3 kj^ 
should not introduce sizeable effects, since the values of the variance in this interval is about 
one order of magnitude smaller than that of the most unstable mode. 



B. Details of the interaction 

In the present paper we adopt the same schematic Skyrme-like effective interaction as in 
Ref. |22|: 

A,,j{k) = A{k) + Sijik) . (5) 
For the symmetric term A{k) we use the finite-range effective interaction introduced in Ref. 
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A{k) = {aI- + {a + l)-^Q'')e-^'''^' , (6) 

Qeq Qeq 

where q = + ^2 is the uniform mean value of the total density of nucleons, Qeq = 0.16 fm~^ 
is the density of symmetric nuclear matter at saturation and 

A = -356.8 MeV, B = 303.9 MeV, a= - . 

6 

The width of the Gaussian in Eq. ([6]) has been chosen in order to reproduce the surface 
energy term as prescribed in Ref. [29]. 

The isospin-dependent part, Sij{k), contains three different terms 

S.Ak) = §^ + r.r.DP + ^-±^Vc{k)5,, , (7) 

with Ti = 1 and T2 = —1. Here Sgym represents the potential part of the symmetry-energy 
density. For the coefficient of the isovector surface term we use the value D = AO MeV fm^ 
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30| 

Si 



Moreover, we include the Coulomb interaction Vc{k) according to the approach of Ref. 
A mean-field exchange contribution 

is also added to the bare Coulomb force. 

In order to stress the effects of the asymmetry of the nuclear medium, we will present 
results obtained with two different parametrizations of the symmetry energy: one with a 
stronger density dependence ( "superstiff " asymmetry term) and the other one with a weaker 
density dependence ("soft" asymmetry term). In both cases the density dependence of the 
potential part of the symmetry-energy density can be expressed by 



SsymiQi, Q2) = S{g){g2 - QiY 



with 



where d = 19 MeV 



S{g) 



2d 



QlQeq 

for the "superstiff" case, and 



S{q) = di — d2Q ■ 



33 , for the "soft" case. 



(8) 
(9) 

(10) 
It should 



where di = 240.9 MeV fm^ and ^2 = 819.1 MeV fm^ 
be noticed that S{g) is nothing but the potential part of the symmetry-energy coefficient 
divided by g, S{g) = C^y^^{g)/g. The inclusion of the Coulomb interaction gives rise to an 
overall decrease of the growth rate of density fluctuations with a corresponding contraction 
of the instability region in the [g, T) phase diagram 3ll, |34 



C. Isospin effects 



Proton and neutron densities oscillate in phase and out of phase, respectively in 
the isoscalarlike fluctuations and in the isovectorlike fluctuations, although with differ- 
ent amplitudes in general. The ratio between amplitudes is given by erf -^{k) / af 2{k) = 

i{k) / (T2 2{k)) , with + for the isoscalarlike case and with — for isovectorlike case. The 
above relation follows from the relevant property 

det\Cl^f{k)\=0, 



(11) 



of the residues. Hence, for asymmetric matter, even in unstable isoscalarlike oscillations that 
lead to phase separation, protons and neutrons move with different amplitude. In particular 
one observes that the ratio between proton and neutron fluctuations is larger than the Z/N 
ratio of the original matter, leading to a more symmetric liquid phase, the so-called isospin 
distillation effect (lol . 11 . 12, 13, IJ, l^- The ratio between the proton and neutron density 
fluctuations is mostly determined by the isoscalarlike fluctuations. This is simply related 
to the fact that during the spinodal decomposition, the isoscalarlike fluctuations are much 
larger than the isovectorlike ones. Thus, this ratio, for a given value of k, can be put as: 



Sgiik) 



Ctik) 



(12) 
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where the ratio of the residues at the imaginary pole, iTk, of Dij{k,u) is given by 



(fc) _ D';'\k,tT,)il- D'^\k,tT,)A2,2{k)) 
CUk) Df\k,iVk){l-Df\k,oj)AiM) 



(13) 



The non-interacting particle-hole propagator, in the semiclassical approximation, is ex- 
pressed by [23] 

z,<«.,Mr,,.-|i + i,„^F(,,,)^. 

where the effective chemical potential fii of neutrons or protons is measured with respect to 
the uniform mean field f/j(f?i, ^2) of the unperturbed initial state, and F{(3jli) is the function 



with j3 = 1/T being the inverse temperature. For the physical parameters considered in 
the present paper the second term in the expression of Df\k,iVk) can be neglected at a 
satisfying approximation, then Eq. ( |T3l) becomes 



(14) 



where / = /(^, Q?)j with Q-i = Q2 — Qi, is the sum of the kinetic and entropy terms of the free- 
energy density, i.e. / represents the free-energy density of a non-interacting two-component 
Fermi gas with effective chemical potentials fti. 

For a qualitative analysis we can neglect the Coulomb interaction. For small asymmetry 
values, we can limit ourselves to consider only the first order term in the expansion of Eq. 
f fT^ in powers of the asymmetry a = g^/ g. In this case the ratio of residues is given by 

CIW ^ . , 1 f.d d'f{g,gs) . ^^dS{g). 

where B{k) represents the sum of the part of interaction of Eq. ([5]) common to both the 
species of nucleons, and of the second derivatives of the free-energy density f{g, g^): 



B{k) = A{k) + 2S{g) + De+ 2 [^3=0 + 



dgl ''^=° dg^ 



Equation (fT5|) shows that for the isospin distillation effect a crucial role is played by the 
derivative of the coefficient of the symmetry energy with respect to the total density. 

With the parameters of the interactions considered here, we get, for systems moder- 
ately inside the spinodal region, a larger distillation effect with the "superstiff" interaction. 
However, it should be remarked that the derivative of S{g) depends on the density. Actu- 
ally at rather low densities distillation effects are expected to be stronger with the "soft" 
parametrization [ssj. 

The symmetry energy also plays an important role in determining isovector fluctuations, 
6g3. These lead to fluctuations in the (iV — Z) content of a given density domain of mass 
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A = N + Z. For finite values of tlie asymmetry a, it is generally not easy to single out 
from the isospin-dependent interaction, Sij{k), a term which can play a conclusive role in 
determining the widths of isotopic distributions. Only for a = and, in addition, neglecting 
the Coulomb interaction, isovector and isoscalar fluctuations are decoupled. In this case the 
fluctuations Sg^, are only due to the isovector modes, and the residues at the real pole a;™ are 
solely determined by the coefficient S{g) in the expression of the symmetry energy of Eq. (jH]). 
For finite values of a, terms containing derivatives of S{g) contribute to the magnitude of the 
residues C™j{k) as well. Generally speaking, one can expect that the strength of the isovector 
pole increases with the coefficient S{g). Hence, in the limit considered here, uj/T » 1, 
isovector fluctuations becomes larger when increasing the symmetry energy coefficient. It 
should be noticed that also for a 7^ the contributions to 5g^ of the isoscalarlike modes still 
tend to cancel out. 

D. Cluster recognition 

We consider a finite cubic lattice of volume V . We assume the density vanishing at 
the surface of the box. The number of nucleons is then fixed and in the box the liquid 
and vapour phases coexist. Moreover, finite size effects can be taken into account in this 
way. However, it should be remarked that the condition of vanishing density at the border 
of the box imposes some symmetries to the problem. In fact, in this case the number of 
independent oscillations is reduced by 2^. Accordingly, due to the imposed symmetries, 
also fluctuations variances are reduced. This is not so important in the isoscalarlike case, 
since these modes are unstable and fluctuations are amplified anyway. On the contrary, 
isovectorlike fluctuations remain quenched (by a factor 8). However, this does not affect 
our conclusions about the relative comparison of the results obtained with the different 
parametrizations of the symmetry energy. 

In order to extract the spatial density 5pi(r), we perform a Fourier expansion of the 
fluctuations (5pj(k). The Fourier components contain the product of three sine functions: 



with ki = 27C rii/L and x = n^b, y = nyb, z = nj) {rti and n^^y^z take positive integer values), 
where L and h are the lengths of the sides of the box and of the primitive cell respectively. 
The coefficients of the functions 0k (i") in the Fourier expansion are linear combinations of 
the isoscalarlike and isovectorlike Fourier coefficients 5^j(k, t), given by stochastic processes 
with the probability distribution of Eq. ([2]) and the ratio of proton to neutron amplitudes 
given by (jl^Ak) / a^^ik) ■ 

The size of the box will be chosen in order that the box contains a number of nucleons 
of the same order of magnitude as in actual heavy ion collision experiments. While the 
adopted size of the primitive cell, Vceii, is such that at least a nucleon of both the species 
is present in the cell on the average. Clusters of nucleons are associated with large density 
domains and will be formed by means of a coalescence algorithm, solely based on the density 
of neighbouring cells. Adjacent cells with a value of the density above (liquid) or below 
(vapour) the average density are collected together. Then, we obtain domains of higher 
density surrounded by domains of lower density. In this way we can investigate separately the 
properties of the liquid phase and of the vapour phase. Moreover, we can build observables 
starting from event by event cluster distributions. 




(16) 
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FIG. 1: (Color online) Charge distributions for clusters belonging to the high density phase for two 
values of the global asymmetry cto, calculated with the "soft" asymmetry term, circles (ao = 0.1) 
and triangles (ao = 0.2), and the "superstiff" asymmetry term, squares (qq = 0.1) and crosses 
(ao = 0.2). 

10"' b' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' H 



10' 



X 
CO 

E 



N 
> 



10 







10 20 30 40 50 



max 



FIG. 2: (Color online) Distributions of the heaviest cluster for two values of the global asymmetry 
uq. Symbols as in Fig. [TJ 



III. RESULTS 



In this section we present and discuss statistical distributions of the clusters obtained 
using the coalescence recipe outhned above. In particular we focus on the isospin content of 
the clusters coming from both the liquid phase and the vapour phase. The value chosen for 
the average density g = OAgeq and the temperature T = 4.5 MeV are in the range expected 
for the multifragmentation process jo], T^]- For the time that the system spends in the 
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FIG. 3: (Color online) Average ratio N jZ as a function of Z for clusters belonging to the high 
density (liquid) phase, calculated with the "soft" asymmetry term (circles and triangles) and the 
"superstiff" asymmetry term (squares and crosses). Bottom: global asymmetry ao = 0.1. Top: 
ao = 0.2. Solid horizontal line: global ratio N jZ for oq = 0.1. Dashed horizontal line: global ratio 
NjZ for ao = 0.2. 



instability region, we have chosen a value of t = 80fm/c. This value is compatible with that 



obtained within the stochastic mean-field approach of Ref. j35|. Moreover, in this short 
time interval the growth of fluctuations is still limited so that a linear approximation can be 
considered as reasonable. For the side of the cubic lattice containing the nuclear system, we 
have adopted the value L = 21 fm, while the side of the primitive cell has a length of 3fm. 
With the chosen value of density the box contains ~ 370 nucleons. 

We have performed calculations for the two parametrizations of the symmetry energy 
introduced above and for two values of the global asymmetry, ao = 0.1, 0.2. We have run 
4 • 10^ events for each case. With the chosen values of the parameters the average numbers 
of nucleons in the liquid phase and in the vapour phase are ~ 270 and ^ 100 respectively. 



A. Charge distributions 

In order to assess the validity of our approach, we have evaluated the distribution of 
the cluster yield. In Fig. [1] we present the charge distribution of clusters belonging to the 
liquid phase, calculated for two values of the global asymmetry (ao = 0.1,0.2) and for the 
two parametrizations of the symmetry energy considered. One can see that the four curves 
displayed in the Figure are rather similar. This is expected from the fact that isoscalarlike 
modes and related variances, that determine the density growth and the appearance of 
fragments, do not much depend on the behaviour of the symmetry energy, nor on the 
global asymmetry of the matter considered. Moreover, the obtained charge distribution is 
similar to the results of full stochastic mean-field simulations [q, i3] of heavy ion collisions at 
~ 30MeV/A, where similar conditions of temperature and density inside the spinodal region 
are encountered. This indicates that fragment size is essentially determined by the properties 
of the most unstable normal modes, as derived in the linear approximation. The beating of 
these modes leads to the rather wide charge distribution. The charge distribution of Fig. [1] 
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FIG. 4: (Color online) The same as in Fig. [3] but for clusters belonging to the low density (vapour) 
phase. 



compares rather well also with the experimental results of Ref. j^, where fragmentation of 
systems of similar size is investigated. 

Another observable of experimental interest is the distribution of the heaviest cluster 
obtained in each event. The study of this observable requires, of course, an event by event 
analysis. Results are shown in Fig. [2l Also in this case we observe that the four curves 
corresponding to different charge asymmetries and to different parametrizations of the sym- 
metry energy look quite similar. Moreover, our results are similar to the predictions of full 
simulations and to experimental data j^. 



B. Isotopic properties 

From Figs. [3] and H] we can appreciate the isospin distillation process occurring in the 
spinodal decomposition of asymmetric nuclear systems. The average ratio N jZ is plotted 
as a function of Z for the correlated density domains of the liquid phase (see Fig. [3]) and of 
the vapour phase (see Fig. Hj). One can nicely see that the ratio N jZ is smaller than the 
initial global value for the liquid, while the opposite holds for the vapour phase. We notice 
that this effect is present also for the two coexisting phases at thermodynamical equilibrium. 



as shown in Refs. |14j . |26j. Moreover, we observe a slight decrease of the ratio N jZ with 
Z . This trend is essentially related to the fact that lighter clusters originate from lower 
density domains where, as discussed in the following (see Subsection C), isospin distillation 
is less effective. We notice that this behaviour, which is in agreement with the results of 



previous calculations based on dynamical models Il4, l36| , is not observed in statistical model 



calculations, that exhibit an opposite trend [37|, |38j. Hence, this property of the ratio NjZ 



with respect to the cluster charge would allow to disantangle between the predictions of 
different models. 

In particular. Figs. [3] and H] show that isospin distillation increases with increasing the 
global asymmetry and is larger when the "superstiff" asymmetry term is used with respect 
to the "soft" case. These features are explained by Eq. (fT5|) . Indeed, the induced asymmetry 
is proportional to the global asymmetry and, at the considered density, the derivative 
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FIG. 5: (Color online) Average distribution of asymmetry a for clusters belonging to the high 
density (liquid) phase, calculated with the "soft" asymmetry term (circles and triangles) and the 
"superstiff" asymmetry term (squares and crosses). Panel (a): global asymmetry ao = 0.1. Panel 
(b): ao = 0.2. 



dS{g)/dg has nearly the same magnitude, but opposite sign for the two used parametriza- 
tions of the symmetry energy (positive for the "superstiff" case and negative for the "soft" 
case) . 

Now we turn our attention to the probabihty distribution of the asymmetry parameter 
for the domains belonging to the two phases. In Figs. and E] we report the probability of 
finding a cluster with the asymmetry a = [N — Z)/ [N + Z) for the liquid and vapour phases 
respectively. The distribution is averaged over intermediate mass clusters (5 < + Z < 50) 
for the high density phase and over light clusters {3 < N + Z < 15) for the low density 
phase. The shift of the maximum toward the left (right) for the liquid (vapour) phase in 
the "superstiff" case with respect to the "soft" case is due to the more effective distillation 
power of the "superstiff" interaction, as already seen above. In addition we observe that 
distributions are wider in the "soft" case. This effect can be related to the fact that at the 
considered density the coefficient S{g) of Eq. (jS]) is larger in the "soft" case with respect 
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FIG. 6: (Color online) The same as in Fig. [5] but for clusters belonging to the low density (vapour) 
phase. 



to the "superstiff" case, leading to a larger value of the residues in Eq. (jll). Moreover, 
the distributions are generally broader for the vapour phase. In the following, we will 
show that this feature can be easily explained on the basis of the relative contributions of 
isoscalarlike and isovectorlike oscillations, at least within our linear treatment of fluctuations. 
Moreover, we observe that the asymmetry distributions become narrower with increasing 



the global asymmetry for both the two phases. This was observed also in Ref. [26|], within 
a thermodynamical study of spinodal decomposition. A final remark on this topic is in 
order. As already seen for the distillation effect, the behavior of the width of the isotopic 
distributions, in connection with the parametrization adopted for the symmetry energy, 
depends on the density value of the matter undergoing spinodal decomposition. With the 
present approach the widths, though slightly, result to be smaller in the "superstiff" case 



than in the "soft" case, whereas in Ref. |22| we have observed an opposite tendency. In 



fact, with the choice of parameters of Ref. 2^ nuclear matter explored a region of higher 
instability (lower density) with respect to the present calculations. Then, the ratio between 
the amplitudes of the isoscalarlike fluctuations (that also contribute to the isotopic variance, 
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due to the dependence of the distillation effect, Sg2/Sgi, on the wave number k) and of the 
isovectorlike fluctuations was quite larger than here. Furthermore, the weighting procedure 
of the different fluctuation modes in Ref. 2^ could likely underestimate the isovectorlike 
fluctuations. 



C. Density dependence of distillation effect 

Now we show that the amplitude of the isospin distillation effect, as well as the variances 
of the isotopic distributions, is related to the value of the density reached in the different 
domains associated with nucleon clusters. 

Due to the beating of the several unstable modes, clusters with a given charge are not 
associated with a fixed density domain, but may come from domains of different volume 
having different density, i.e. the cluster density g may fluctuate. We will show that this 
effect influences the cluster isotopic properties. 

As far as the distillation effect is concerned, we notice that the N/ Z ratio of a cluster can 
be written as: 

N/Z={gl + Sg2)/{g1 + 6g,), 

where ^2 are the neutron and proton densities of the initial matter, while 5g2 and 

5gi are the corresponding fluctuation values in the domain considered. After some algebra 
one obtains: 

^j^ly, [{NlZ),-5g2/5g,]-5g,/g\ 
N/Z - {N/Z), 

where {N/Z)q = g^/ gi- We notice that, due to the distillation effect, Sg2/Sgi < {N/Z)q. 
From the expression given above, one can see that the ratio N/Z of a given cluster decreases 
when the density g of the domain increases, i.e. 5gi (and 6g2) are larger. In fact, larger den- 
sities are associated with larger isoscalarlike fluctuations and thus with a larger distillation. 
On the other hand, clusters in the vapour phase, having density smaller than the initial one 
{6g < 0), show a larger value of the ratio N/Z. Moreover, the deviation of the ratio N/Z 
with respect to the initial value {N/Z)o increases as far as the ratio 5g2/Sgi gets smaller, 
i.e. for parametrizations of the symmetry energy that lead to larger distillation effect. In 
Fig. [7] we show the N/Z ratio, averaged over all charg GS, clS db function of the cluster den- 
sity. The decreasing trend is clear, especially for the more neutron-rich system, where one 
can also better appreciate the different effect of the two equations of state employed in the 
calculations. 

For the variance of the isotopic distributions, we also expect a decreasing behaviour with 
density. In fact isovectorlike fluctuations, that are mainly responsible of the isotopic variance, 
in our linear approach are implemented according to the properties (density, asymmetry, 
temperature) of the initial matter. Then, nucleon clusters with larger density should have 
narrower isotopic distributions due to the larger contribution of the isoscalarlike modes with 
respect to isovectorlike fluctuations. 

The probability distribution for the isovector density fluctuations 6p3 in a given domain 
of volume V can be approximatively expressed as: 

P^exp{-6pl/2al), (17) 

where the variance af^ can be evaluated starting from the value of crfjik) of Eq. (jlj) and 
is inversely proportional to the volume V, af^ = g{g,T)/V, where g{g,T) is a function 
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FIG. 7: (Color online) Average ratio N/Z as a function of the cluster density for two values of the 
global asymmetry ao, calculated with the "soft" asymmetry term, circles (oq = 0.1) and triangles 
(ao = 0.2), and the "superstiff" asymmetry term, squares (ao = 0.1) and crosses (ao = 0.2) 



that depends on the initial conditions of the matter undergoing spinodal decomposition and 
on the parameters on the interaction adopted [3]. Then, for a cluster of volume V, the 
distribution P{N — Z) can be written as: 

P{N -Z)^exp(-[N-Z-{No- Zo)]' /[2Vg{g, T)]). (18) 

Hence, the variance of the isotopic distribution should be proportional to V. This is con- 
firmed by our calculations, as shown in Fig. [HI where we plot the variance of the {N — Z) 
distribution, averaged over all clusters with volume V = NyVceiu ^-s a function of Ny. From 
Fig. [S]one can also see that, especially in the case of the "superstiff" interaction, the vari- 
ance is smaller for the neutron richer system (compare squares and crosses). Moreover, for 
a given volume, variances are larger in the "soft" case. 

If we consider now clusters with a given mass A, on can easily realize that the isotopic 
distributions will depend on the cluster density, being broader for clusters with smaller 
density (i.e. larger volume). Hence clusters that originate from larger density domains, are 
expected to have an average N/Z ratio smaller than the value corresponding to the initial 
matter and, at the same time, are expected to show narrower isotopic distributions. 



IV. CONCLUSIONS 

In this paper we have considered a piece of asymmetric nuclear matter contained in a 
cubic box with impenetrable walls. The values of density and temperature of the system 
are inside the spinodal instability region of the phase diagram. We have mainly focused our 
attention on the isotopic distributions of clusters coming from the break-up of the system. 

In our approach the density fluctuations are considered within a linear approximation, 
then our results refer to nuclear matter being not deeply inside the unstable spinodal region. 
Such physical situation can occur in central collisions of heavy ions at moderate energy 
(~ ?){)MeV / A). Moreover, we remark that our investigation concerns the distributions just 
after the early break-up of the nuclear system. 
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FIG. 8: (Color online) Average variance of the (A^ — Z) distribution for clusters of volume V = 
NyVceii as a function of the number of the elementary cells Ny calculated for two values of the 
global asymmetry oq. Symbols as in Fig. [71 



It has been found that, at the considered values of density and asymmmetry, the global 
features of the disassembly, the isospin distillation effect included, are essentially determined 
by the unstable isoscalarlike fluctuations of the density, while the widths of the isotopic 
distributions are mostly affected by the isovectorlike oscillations. 

Within the model developed it is possible to recognize clusters originating from domains 
having different density values. Discerning clusters with density lower or higher than the 
initial density, we have observed a clear occurrence of the isospin distillation effect: low- 
density clusters are more neutron-rich than the initial system, correspondingly the opposite 
happens for higher density clusters. As a general trend the average value of the N/Z ratio 
decreases with increasing the cluster density. This effect results to be enhanced for more 
charge-asymmetric systems. Moreover, it has been found that the isotopic distributions 
of clusters belonging to the low density phase are broader than those of clusters belonging 
to the high density phase. Hence larger distillation effects are associated with narrower 
isotopic distributions, at least within our linear treatment of fluctuations. This feature, 
if it survives the secondary decay of primary fragments, could be qualitatively checked in 
experiments by comparing the isotopic distributions of light clusters, that should emerge 
from the low-density phase, and intermediate mass fragments, most likely associated with 
higher density regions. However, at the present level of investigation, a careful comparison 
with experiments would require a reconstruction of the properties of the primary fragments, 
based on actual experimental data [39] . 

In order to stress the sensitivity of the results to the isovector part of the nuclear inter- 
action, we have presented calculations obtained with two different parametrizations of the 
symmetry energy. For the values of charge-asymmetry considered, we have seen that isospin 
effects are essentially related to the values of the coefficient of the symmetry energy and of 
its derivative with respect to the density. In particular, the former quantity mainly deter- 
mines the width of the isotopic distributions, while the latter quantity plays a crucial role 
for determining the strength of the isospin distillation. However, it should be remarked that 
the symmetry energy properties depend on the density reached by the system during the 
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decomposition process. Then, the effectiveness of a given parametrization of the symmetry 
energy for the isospin distillation depends on the actual path followed by the nuclear system 
in the density-temperature phase diagram. 
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